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Abstract. 

We present the status of the joint search for gravitational waves from 
inspiraling neutron star binaries in the LIGO Science Run 2 and TAMA300 Data 
Taking Run 8 data, which was taken from February 14 to April 14, 2003, by the 
LIGO and TAMA collaborations. In this paper we discuss what has been learned 
from an analysis of a subset of the data sample reserved as a "playground" . We 
determine the coincidence conditions for parameters such as the coalescence time 
and chirp mass by injecting simulated Galactic binary neutron star signals into 
the data stream. We select coincidence conditions so as to maximize our efficiency 
of detecting simulated signals. We obtain an efficiency for our coincident search of 
78%, and show that we are missing primarily very distant signals for TAMA300. 
We perform a time slide analysis to estimate the background due to accidental 
coincidence of noise triggers. We find that the background triggers have a very 
different character from the triggers of simulated signals. 



PACS numbers: 95.85.Sz, 04.80.Nn, 07.05Kf, 95.55Ym 

1. Introduction 

In the past several years, there has been substantial progress in gravitational 
wave detection experiments by the ground-based laser interferometers, LIGO PP, 
TAMA300 Hj, GEO600 3 and VIRGO @|. 

The LIGO and TAMA collaborations are conducting joint searches for 
gravitational waves in the LIGO Science Run 2 (S2) and TAMA300 Data Taking Run 
8 (DT8) data, which was taken from February 14 to April 14, 2003. Three classes of 
gravitational wave transients are being sought: unmodeled gravitational wave bursts 
from the gamma ray burst event GRB 030329 5 ; unmodeled gravitational wave bursts 
without an electromagnetic trigger 6 ; inspiral signals from Galactic binary neutron 
star (BNS) systems. 

In this paper, we present the status of the joint search for gravitational waves from 
inspiraling neutron star binaries by the LIGO and TAMA collaborations. We focus on 
what has been learned from a subset (10%) of the full coincident data set, chosen to 
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be representative of the full data set in term of detector contribution and noise, and 
reserved for tuning of the analysis pipeline. We refer to this as the playground data. 

We use the observable effective distance to characterize the sensitivity of the 
instruments to binary neutron star inspirals. This is defined as the distance at which 
an inspiral of 1.4-1.4M© neutron stars, in the optimal direction and orientation with 
respect to each detector, would produce a signal to noise ratio (SNR) of 8. The effective 
distance is always greater than or equal to the actual distance, and on average is 2.3 
times as large as the distance. In Table we show the observable effective distance 
of the three LIGO instruments (HI, H2 and LI) and the TAMA300 instrument (Tl) 
during S2 and DT8. Since the LIGO and TAMA300 detectors were sensitive to the 
majority of the Milky Way, it is possible for us to detect an inspiral signal from our 
Galaxy. If we do not make a detection, then we can place an upper limit on the binary 
neutron star inspiral rate in the Milky Way. Another motivation for this work is to 
gain experience in performing coincident analyses between different groups, and to 
establish the analysis method. 



Detector 


Observable effective distance 


Comments 


LI 


- 2000 kpc 


Reach M31, M32, M33, MHO 


HI 


- 900 kpc 


Barely reach M31 etc 


H2 


- 600 kpc 


Covered the Milky Way. 


Tl 


- 50 kpc 


Covered most of the Milky Way. 



Table 1. The observable effective distance, during the S2 and DT8, at which 
an inspiral of 1.4-1.4 Mq neutron stars, in the optimal direction and orientation 
with respect to each detector, would produce the SNR of 8. 

Each group analyzes its own data, generating a "trigger" when an inspiral 
candidate appears in the data stream. The groups then exchange triggers and search 
for coincidences. Both LIGO and TAMA300 have performed independent analyses on 
the S2 /DT8 data set and they have been described in detail elsewhere |S] . In this 
paper, we will focus on aspects of the coincident search. Using the playground data, 
we define a set of coincidence requirements between triggers from the two detectors. 
Before we perform the coincident analysis using these coincidence conditions, they are 
tested by performing software injections of simulated inspiral signals into the data and 
by checking the detection efficiency. We verify that the injection efficiency is consistent 
with the observable effective distances stated in Tabled Finally, using a technique 
of artificially shifting the time series with respect to one another, we estimate the 
coincident triggers produced accidentally by the instrumental noise. 

This paper is organized as follows. In section |3 we describe the LIGO and 
TAMA300 detectors and the data sets to be analyzed. In section|31 we briefly describe 
the single instrument playground analysis. In section 0] we discuss the method of 
coincident analysis using the results of the single-detector searches. In sectional we 
estimate the background of coincident triggers using a time slide analysis and test the 
efficiency of our search using simulations. Section [S] summarizes the paper. 

2. Detectors and data sets 

The LIGO network of detectors consists of a 4km interferometer "Ll"in Livingston, LA 
and a 4km "HI" and a 2km "H2" interferometer which share a common vacuum system 
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in Hanford, WA. TAMA300 is a 300m interferometer "Tl" in Mitaka, Tokyo. Basic 
information on the position, orientation of these detectors and detailed descriptions 
of their operation can be found in 12 ■ 

The data analyzed in this search was taken during LIGO S2, TAMA300 DT8, 
between 16:00 UTC 14 February 2003 and 16:00 UTC 14 April 2003. We only analyze 
data from the periods when both LIGO and TAMA300 interferometers were operating. 
Furthermore, we restrict to times when only one of the LIGO sites was operational, 
thus sharing no common observational time with the LIGO only search discussed in 

Therefore, we have four independent data sets to analyze: the data set during 
which neither HI nor H2 were operating, the nHl-nH2-Ll-Tl coincident data set (here 
"n" stands for "not operating" ) ; and three data sets when one or both of the Hanford 
detectors were operational but LI was not — the Hl-nH2-nLl-Tl coincident data set; 
the nHl-H2-nLl-Tl coincident data set; and the Hl-H2-nLl-Tl coincident data set. 
The total length of data was 650 hours. This comprised 334 hours of H1-H2-T1 data, 
212 hours of Hl-Tl data, 68 hours of H2-T1 data and 36 hours of Ll-Tl data. To 
avoid any bias from tuning our pipeline using the same data from which we derive 
our upper limits, the tuning of analysis parameters was done without examining the 
full coincident trigger sets. Instead, preliminary tuning was done on the playground 
data, which will not be used placing the upper limit. In this analysis, the length of 
playground data was 64 hours. 

3. Overview of the single interferometer analysis 

Our analysis methodology is similar, though not identical, to that used in the LIGO 
analysis [5] [U] and the TAMA300-LISM coincident analysis JU| ■ We analyze the data 
from each detector independently, identifying candidate inspiral events in the data 
streams which we refer to as "triggers" . We then search for coincidences within the 
lists of triggers from each interferometer. Here, we will only briefly review the main 
steps in the single instrument analyses. The search methods employed by TAMA300 
are discussed in greater detail in 0E3 an d by LIGO in piTT]. 

For each instrument we produce a bank of inspiral templates with the 2PN 
expansion for LIGO [T2\ and the 2.5PN expansion for TAMA300 [Hj, distinguished by 
the masses (mi and m-i) of the two neutron stars in the binary. We then filtered the 
data stream through this bank of templates and record a trigger whenever the SNR 
exceeds 7. In order to minimise the number of spurious triggers, we also perform a 
test of waveform consistency (the \ 2 test EH)- Additionally, during times when both 
Hanford detectors are operational, we check for coincidence between HI and H2. In 
the playground data, there were no such triggers. However, if a trigger is seen in HI 
and its amplitude, or estimated effective distance, is such that we would not expect 
to see it in H2, we keep this HI trigger even with no coincident trigger in H2 (see [H] 
for further details). 

The triggers from each instrument are characterized by the time of coalescence t c , 
the chirp mass M. = Mr] 3 / 5 , the symmetric mass ratio (dimensionless) rj = ■m-ira^jM' 1 , 
(M = mi + rri2 is the total mass of the binary system), the signal to noise ratio p, 
value of the \ 2 test and effective distance D c s . The effective distance of a waveform is 
the distance for which a binary neutron star system would produce the waveform if it 
were optimally oriented. The lists of triggers are then exchanged between the LIGO 
and TAMA analysis groups and coincident triggers are searched for. 
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4. Tuning of the coincident analysis 

True gravitational wave events, if sufficiently loud, will appear in the trigger lists 
of both the LIGO and TAMA300 instruments. However, they will have somewhat 
different values of parameters such as coalescence time and chirp mass. This is due to 
the detectors' noise, the difference in the detectors' locations and arm orientations, and 
the discrete sampling of the mass parameter space by the template bank. To evaluate 
the accuracy with which we can determine various parameters, we perform Monte 
Carlo simulations of injected signals. Our sample population for simulations consisted 
of Galactic neutron star binaries with masses in the range I.QMq < mi, m 2 < 3.0M Q . 
For the simulations, the spatial and mass distributions were taken from jf 5} and |16|. 
We set coincidence windows to maximize our detection efficiency while reducing the 
number of chance coincidences. In the LIGO-TAMA300 search we require coincidence 
in both time and chirp mass. 

4-1- Timing accuracy 

We require the coalescence time t c of a signal to be consistent between the LIGO and 
TAMA300 detectors. The allowed time difference is determined by our accuracy in 
determining t c and the light travel time between the sites. To determine our timing 
accuracy, we inject simulated gravitational wave signals into the data streams from 
each detector, and re-analyze the data. We injected 660 signals in TAMA300 data 
and 455 signals in LIGO data. We compared the injected time and the detected time 
which was obtained from our search pipeline. We define the detected time to be the 
time of the loudest trigger within 3 ms of the coalescence time of the injection. We 
show the accuracy of coalescence time; detected - injected coalescence time in Figs. ^ 
and [21 We find that nearly all injections into both LIGO and TAMA300 data have 
triggers (p > 7) recorded within ±1.0 ms of the injection time. 



Histogram of End Time Accuracy 




Detected - Injected E nd Time (msec) 

Figure 1. Accuracy of the reconstructed coalescence time t c for simulated 
injections into LIGO data. All of the 455 injected signals were seen (p > 7) 
in the LIGO data, and the coalescence time of the majority of them was recorded 
within ±1.0 ms of the injection coalescence time. 



The distance between the TAMA300 site and the LIGO Hanford site is 7487 
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Histogram of End time Accuracy 
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Detected - Injected End Time (msec) 

Figure 2. Accuracy of t c for simulated injections into TAMA300 data. Out of 
the 660 injections into TAMA300 data, 516 produced triggers (p > 7). For the 
majority of these injections, the coalescence time was recorded within ±1.0 ms of 
the injection time. The other simulated signals injected into TAMA300 data did 
not produce triggers because they were too weak to be detected by TAMA300. 



km. The distance between the TAMA300 site and the LIGO Livingston site is 9703 
km. Because gravitational waves travel at the speed of light, the maximum delay of 
the arrival time of gravitational wave signals are A^istH-T = 25.0ms for TAMA300- 
Hanford and A£<j ist L_T = 32.4ms for TAMA300-Livingston. Based on these facts, we 
define the allowed difference in coalescence time At w by the maximum time delay of 
signal plus the uncertainty in determining the coalescence time. The values are 
At w h-t = 27ms for Hanford-TAMA300 and 

At w l-t = 35ms for Livingston-TAMA300 . (1) 
Therefore, if the coalescence times, £ Ci ti for TAMA300 and t c ,i (where i = H or L) 
for LIGO, of a trigger satisfy, 

\tc,Tl — t Cy i\ < Atw i-T, (2) 

the trigger is recorded as a candidate event. 



4-2. Chirp mass accuracy 

The accuracy with which we can determine the chirp mass is determined in the same 
way as for the coalescence time; namely by adding simulated signals into the data of 
the instruments. As with the end time, we take the detected chirp mass to be the chirp 
mass of the loudest trigger obtained within 3 ms of the end of the injection. In Figs. |3] 
and^J we find that the chirp mass of an injected signal is recorded with an accuracy 
of 0.01M© in the LIGO instruments and 0.05M Q in TAMA300. We define the allowed 
difference in chirp mass of coincident triggers between LIGO and TAMA300 to be: 

AM W = 0.05M Q . (3) 

If the chirp mass, Mti for TAMA300 and Mi (where i = H or L) for LIGO, of a 
trigger satisfies, 

\Mn-Mi\<AM w , (4) 
the trigger is recorded as a candidate event. 
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Figure 3. Accuracy of chirp mass M for simulated injections into LIGO data. 
The majority of the 455 injections into LIGO data have triggers (p > 7) recorded 
with chirp mass within O.OIMq of the chirp mass of the injected signal. 



Histogram of Chirp Mass Accuracy 




Detected - Injected Chirp Mass (M ) 



Figure 4. Accuracy of chirp mass A4 for simulated injections into TAMA300 
data. Out of the 516 found injections, the majority have triggers (p > 7) recorded 
with chirp mass within 0.05Mq of the chirp mass of injected signal. 



5. Background, efficiency and consistency check of the analysis 

Before doing the coincidence search, we study the background events by performing 
a time slide analysis of the data. Additionally, we evaluate the search efficiency by 
performing an injection analysis. We check that the injection efficiency is consistent 
with the sensitivities of the instruments reported in Tabled 

5.1. Background 

Even in the absence of gravitational wave signals, we expect that some noise triggers 
will be coincident by chance. We can estimate the number of such accidental coincident 
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Figure 5. The signal to noise ratios Pligo vs Ptama of the accidental coincident 
triggers, in playground data, using 100 time slides. 



triggers by performing our coincidence analysis on the data with artificial relative time 
shifts added [TJ\. We performed 100 time slides of ±5, 10, 15, 20, ... , 250 sec of the 
LIGO triggers relative to the TAMA300 triggers. These slides are much longer than 
the light travel time between the sites, so that any coincidence cannot be from actual 
gravitational wave. They are also longer than the detector noise auto-correlation 
time, longer than the longest signal template duration and shorter than timescales of 
detectors' non-stationarity, so that each time slide provides an independent estimate 
of the accidental coincident rate. Additionally, the triggers are clustered over a 100 
ms window. The signal to noise ratios p LIGO vs /) TSMA of the accidental coincident 
triggers found in the playground data is shown in Fig EI From these time slides, we 
obtain an estimated false alarm rate of 0.4 triggers per day. 



5.2. Efficiency and consistency check of the analysis 

We performed a set of coherent injections into the data of both TAMA300 and the 
LIGO instruments in order to test our coincidence windows and determine the search 
efficiency. Overall, this consisted of 381 injections in times when TAMA300 and one 
of the LIGO sites were operational. These injections were same as those described in 
Section 01 used for determining the coincidence parameters. 

Among the 381 injections, 297 events produced coincident triggers in the two 
instruments. Of the 84 injections which were not found by the coincidence analysis, 
only one injection was found separately by both LIGO and TAMA300. However, 
in TAMA300, the chirp mass was recovered with an error of 0.07M Q , while LIGO 
recovered it with an error of less than 0.01M Q . Hence, although a trigger was found 
in both instruments, it did not pass the chirp mass coincidence requirement. 

The other 83 missed injection signals were not detected by TAMA300 with p > 7 
and x 2 below threshold. We have checked that the missed injections are actually very 
distant signals for TAMA300. To do this, we plot the effective distance versus GMST 
for the found and missed injections in Fig. [BJ It is clear that the majority of missed 
injections are too distant for TAMA300 to detect. We also plot the efficiency versus 
effective distance in Fig[7| We find that the efficiency is very close to 100% within 
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Plot of found and missed injections against TAMA effective distance 
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Figure 6. Injected distance versus GMST for Galactic BNS inspiral signals 
injected into the TAMA300 playground data stream. This figure shows injected 
signals which were found by the LIGO-TAMA300 inspiral pipeline (x), and those 
which were missed (+). The majority of missed injections are at distances greater 
than the range of the TAMA300 instrument during DT8. 



Plot of LIGO - TAMA efficiency against TAMA effective distance 
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Figure 7. LIGO-TAMA300 coincident detection efficiency for Galactic BNS 
inspirals versus effective distance of the signal at the TAMA300 detector. 



50 kpc and then drops off steeply. This is consistent with the observable distance of 
TAMA300 given in Table □ 

In Fig. we show the scatter plot (/9 L igo> Ptama) with both the time slide and 
injection triggers from the playground analysis. The plot shows clearly that the 
background triggers are well separated from injections of Galactic BNS inspiral signals 
making it easy to distinguish the latter from the former. Therefore, it would be possible 
for us to distinguish a Galactic inspiral signal from the background noise. 
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Figure 8. The signal to noise ratios of inspiral triggers from LIGO versus 
TAMA300, including the time slide triggers (•) and the playground injection 
triggers (x). 



6. Summary 

We presented the status of the joint search for gravitational waves from inspiraling 
neutron star binaries by the LIGO and TAMA collaborations. We focused primarily 
on what has been learned from playground data. 

We discussed a trigger based method of performing the coincidence analysis in a 
multi detector search for gravitational waves from inspiraling neutron star binaries. 
We determined efficient coincidence conditions on the coalescence time and chirp 
mass by analyzing Galactic binary neutron star signals injected into the data streams 
of the LIGO and TAMA300 detectors. The choice of these coincidence parameters 
was validated by performing a set of coherent simulated injections into both LIGO 
and TAMA300 playground data. Using these injections, we found the detection 
efficiency to Galactic BNS inspirals for our coincident search to be 78%. This value 
was consistent with expectations based on the sensitivities of the instruments during 
S2/DT8, derived from their typical noise spectra. Additionally, we performed a time 
slide analysis of the playground data. This allowed us to estimate our expected 
background due to chance coincidences. Furthermore, we saw that the character of 
the injected signals was very different from that of the background triggers. 

The coincident analysis of all of the LIGO S2-TAMA300 DT8 coincident data 
will be performed using the methods discussed in this paper. Complete results of the 
analysis which will include any detection candidates as well as an upper limit on the 
event rate obtained from all of the LIGO S2-TAMA300 DT8 coincident data will be 
reported in a subsequent publication. 
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